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Abstract 

Phonons are universal language for solid state theory. Effective phonon theory (self-consistent 
harmonic approximation) have been extensively used to access the weakly nonlinear effect in solid 
materials, which are standard content in textbook. In this manuscript, we test the effective phonon 
theory in one-dimensional chains with asymmetric interactions. It is found that simulation results 
agree well with theoretical predictions for symmetric cases, while significantly deviate from theo- 
retical predictions for asymmetric cases. Our results imply the asymmetric interaction can provide 
stronger phonon-phonon interaction, which lead to break down of the effective phonon theory. 
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Nonlinearity play important roles in solid materials, which responds for thermodynamic 
and transport properties such as thermal expansion and thermal conductivity. Effective 
phonon theory (EPT) have been commonly used to assess the weakly nonlinear effect on 
phonons (normal modes in classical physics), whose basic ideal is to find a effective har- 
monic Hamiltonian to approximate the true nonlinear Hamiltonian (in early literature, also 
called self-consistent harmonic approximation) [1]. In the last decade, EPT have been exten- 
sively used to study low dimensional lattice models. These studies focus on one- dimensional 
(ID) momentum-conserved lattices and reported that the numerical results agree well with 
theoretical predictions even in strong nonlinearity limit 2j-l6|. It seemingly suggest EPT 
worked well for ID lattice systems. However, Our investigation here shows things are not 
that simple. We find that simulation results agree well with theoretical predictions for ID 
chains with symmetric nonlinear interactions, while obviously deviate from predictions for 
ones with asymmetric interactions. 

A ID momentum-conserved lattice are generally defined by the dimensionless Hamilto- 
nian 

ff = Ei+%-^-0 (i) 

where pi and Xi are the zth particle momentum and displacement from the equilibrium 
position, respectively. V is the potential between two neighboring particles, which is given 
by 

V(x) = \x 2 + |x 3 + ^ (2) 

The lattice spacing is set to unity. This means the sound velocity c s in the harmonic limit is 
equal to one in our models. The EPT calculate the effective particle force constants through 
self-consistently replacing the harmonic force constants by their thermal averages over all 
possible motions of particles. As a result, an effective harmonic potential approximate the 
nonlinear one (2) by 

V{x) = |* 2 (3) 

where rj is the renormalization factor. To do so, the phonon frequencies are renormalized to 
Uk = V u k by the factor r/ as well, where k is mode number and u is the phonon frequency 
for (2) in harmonic limit, r/ is usually considered as k- independent, which is numerically 
confirmed in FPU-/? models (2, 5|. Many theoretical methods were used to deal with model 
(2) to calculate renormalized factor rj, such as generalized virial theorem 

as 

, Zwanzig-Mori 
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projection |4|], weak turbulence theory [5fl and self-consistent phonon theory |6j|. It have been 
confirmed that these methods are equivalent 4j-[6j. In the present manuscript, we test two 



two kind of expressions for r\ based on the theories mentioned above. One is given by 
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where K and Uh is the total kinetic and harmonic potential energy of model (2). Another 
slightly different expression was given by ^ 



V 



where 5i = Xi — x^\. < • • ■ > in the above two expressions denotes averaging over the Gibbs 
measure. 

We divide the potential fl2j) into two types according to its symmetry with respect to its 
zero point at x — 0. In the symmetric type, we fix a = 0, which recovers the Fermi-Pasta- 
Ulam (3 (FPU-/3) models. In the asymmetric one, we fix = 1 and varied a in the interval 2 
to ensure only one potential minimum at x — 0, which is the Fermi-Pasta-Ulam a/3 (FPU- 
a/3) models. Nonlinear parameter a also governs the degree of the interaction asymmetry. 
By increasing a from zero where the potential is symmetric, one gets increasingly stronger 
asymmetry. 

The effective sound velocity c s is an important outcome from EPT, which can be mea- 
sured in a laboratory. Therefore, it is natural to consider c s as a observable to verify the 
effectiveness of EPT. Obviously, c s equals the renormalization factor rj in our dimensionless 
models. On the other hand, we can employ equilibrium fluctuation correlation method Q 
(for an more detailed version, see [8]) to directly determine c s by simulations. The method 
has also been used to measure the effective sound velocity in FPU-/3 models and the re- 
sults agree well with the theoretical predictions from Eqn. (5) [9]. We conduct numerical 
simulations at constant energy by integrating the equations of motion governed by (2) with 
periodic boundary conditions. We use random initial conditions by assigning velocities from 
a Gaussian distribution and zero relative displacements for each particle, respectively, with 
the two constraints that the total momentum is zero and the total energy is set to be a spec- 
ified constant. The Runge-Kutta algorithm of 7-8th order is adopted with the time step 0.01 
to ensures the conservation of the total system energy and momentum up to high accurate 
in the all run time. The system then evolves for a certain transient time in order to start the 



measurements in well-defined equilibrium states. In order to confirm that the system has 
reached the thermal equilibrium state, we have checked that the value of energy localization 
to be 0(1) and the relative displacement and the velocity distributions of particles agrees 
well with the Gibbs statistics j|3|. 

In Fig. 1, we compare the measured c s as a function of nonlinear parameters with the 
theoretical predictions from Eqs. (4) and (5) where < • • • > is replaced by time average. 
It is clearly seen that the measured results agree well with the predictions of EPT for 
FPU-/3 chains in the entire range of ft under study. In contrast, the systematic deviation 
from the predictions is clearly revealed as a increasing for FPU-a/3 chains. The theoretical 
results obviously underestimate the effective sound velocity of FPU-a/3 chains with larger 
asymmetry. These results imply that EPT does well for nonlinear chins with symmetric 
interactions, but break down for ones with asymmetric interactions. 

For understanding the breakdown of EPT in ID chains with asymmetric interactions, we 
shall perform a more microscopic measurement as follows. We directly extract the phonon 
frequency uj m from the true dynamics governed by (2), then calculate the measured renor- 
malization factor by r\ m = — , where Uh is the corresponding harmonic phonon frequency. 
To this aim, we shall analyze the power spectrum of velocity of a singe particle randomly 
picked up from the chain at equilibrium. The same simulation process as mentioned above 
is used. A FFT algorithm is used to obtain the power spectrum with a time series of total 
time 10 5 to ensure the high frequency resolution. The power spectrum is averaging over at 
least 20 trajectories to suppress the statistical fluctuations. 

Figure. 2 shows the power spectrum of the single-particle velocity time series for FPU- 
(3 model with (3=10 and FPU-a/3 model with a=2. In spite of strong nonlinearity and 
long system size (N = 2048 for both models), it is clearly seen that a few phonon peaks 
with the lowest frequencies is well distinguished from others whose peaks quickly decay into 
noise signals due to stronger the phonon-phonon interactions. It was commonly accepted 
that the long- wavelength phonons dominate the asymptotic behavior of ID momentum- 
conserved lattices. This is the very reason that a long- wavelength approximations was 
extensively adopted in varied theories, such as the self-consistent mode-coupling theory 



and hydrodynamics 
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14j . In addition, some important physical quantity, such as sound 



velocity, is defined in a long-wavelength limit. Therefore, it is fair to focus on the lowest 
frequency phonon peak. The separability of the lowest frequency phonon peak have been 
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FIG. 1: The effective sound velocity is plotted versus nonlinear parameters f3 in FPU-/3 chains (a) 
and a in FPU-a/3 chains (b) at equilibrium temperatures T = 0.5. The open triangles indicate the 
measured value by the equilibrium fluctuation correlation method for N = 2048 7]. The red lines 
indicate the predictions from the Eqn. (4) and the blue lines indicate the predictions from the Eqn. 
(5), where N = 128. 

checked for all parameters under investigation, and hence allow us to accurately measure 
its frequency and renormalization factor. Interestingly, our results intuitively display the 
image of "effective phonons" that was proposed on purely phenomenological basis in varied 
theories. 

Figure. 2 also shows a notable difference between FPU-a/3 and FPU-/3 models for 
linewidths of phonon peaks. The phonon peak for FPU-a/3 model is broader than the one 
for FPU-/3 models by one order of magnitude under the current parameters. That means 
the asymmetric interactions provide stronger phonon-phonon interaction compared to the 
symmetric ones, which lead to shorter phonon lifetime. The results on phonon linewidths 
will be reported in our otherwise manuscript. Here, we focus on the effect of nonlinearity 
on phonon frequency. 
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FIG. 2: The phonon peaks in the low frequency regim for N = 2048and T = 0.5. Blue solid line 
indicate the FPU-/3 with /3 = 10; Red solid line indicate the FPU-a/3 with a = 2 and (3 = 1. For 
comparison, the lowest harmonic frequency is also plotted by black solid line. 

In Fig. 3, we plot the measured renormalization factor i] m as the function of nonlin- 
ear parameters for two representative equilibrium temperatures 0.1 and 0.5, respectively. 
Meanwhile, we calculate the theoretical value of rj via the expressions (4) and (5), in which 
< • • • > was replaced by time average. For FPU-/3 chains, the measured factor i] m agree 
well with theoretical predictions in the entire parameter region being explored at both tem- 
peratures. In sharp contrast to that, the systematic deviation from the prediction is clearly 
seen again for FPU-a/3 chains at larger a at both temperatures. We also check that the rj m 
is independent of systems size N for both FPU-/? and FPU-a/3 models, as shown in Fig. 4. 
This result implies that the measured r] m is independent of the wave numbers of phonons 
for the two models. Therefore, It is reasonable to measure the renormalization factor via 
the single phonon peak with the lowest frequency. The results further reveal that the EPT 
fail to renormalize the phonon frequency in chains with asymmetric interactions. 

Naturally, a following question is if the measured renormalization factor can exactly 
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FIG. 3: The renormalization factor is plotted versus nonlinear parameters (3 in FPU-/3 chains and 
a in FPU-q/3 chains at equilibrium temperatures T = 0.5 and T = 0.1. The open circles indicate 
the measured value T] m . The red lines indicate the predictions from the Eqn. (4) and the blue lines 
indicate the predictions from the Eqn. (5). 

predict the effective sound velocity for the chains under investigation. To check this, we 
compare the effective sound velocity obtained by rj with the measured one from the equilib- 
rium fluctuation correlation method mentioned above. In Fig. 5, it is clearly seen that there 
is an excellent agreement between these two ways to obtain the effective sound velocity in 
both FPU-/3 abnd FPU-a/3 models. These results confirm that one can obtain the correct 
renormalized frequency of phonon even in chains with larger interaction asymmetry where 
the EPT break down. 

The effectiveness of EPT for FPU-/3 chains even in a strongly nonlinear regime has also 
been confirmed by other authors , and be well understood under the frame of weak 

turbulence theory [5]. In these theoretical scenario, the FPU-/3 dynamics with the strong 
nonlinear limit in thermal equilibrium can be effectively described by a system of weakly 
interacting renormalized normal modes. Such effective renormalization results mainly from 
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FIG. 4: The measured renormalization factor is plotted versus system size N at equilibrium tem- 
perature T = 0.5 for FPU-/? chains with (3 = 10 (open circles) and FPU-a/3 chains with a— 2 and 
j3 = l(open squares), where two horizontal lines are responding to their mean values, respectively. 

the trivial resonant wave interactions, i.e., interactions with no momentum exchange. It is 
thus approprite that fro FPU-/3 chains the mean-field approximation is employed as does 
EPT. 

However, such a picture seemingly can not extend to the chains with asymmetric inter- 
particle interactions, as revealed by the results mentioned above. The deviation from the 
calculations in EPT in FPU-a/3 chains suggest that the asymmetric interactions may provide 
stronger nonlinearity than symmetric ones. Mulch-phonon processes beyond the trivial 
resonant wave interactions should be considered, lead to breakdown of EPT. In fact, it was 
also found that asymmetric interactions results in stronger phonon-phonon interactions in 
the recent studies of energy equipartition and heat conduction [isl, in FPU-a/3 chains. 
In the former case, shorter equipartition time was observed for FPU-a/3 chains comparing 
to the FPU-/3 chains, and in the latter case, normal heat conduction was reported, against 
the abnormal heat conduction for FPU-/3 chains. 
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FIG. 5: The measured sound velocity is plotted versus nonlinear parameters /3 for FPU-/3 chains 
(a) and FPU-a/3 chains (b) at the equilibrium temperature T = 0.5. The open circles indicate the 
values obtained by r/ m based on a analysis of the lowest frequency phonon peak for N = 128. The 
solid triangles indicate the values from measured values from the equilibrium fluctuation correlation 
method for N = 2048 [3|. 

In summary, we check the EPT for bath FPU-/3 and FPU-a/3 models through an analysis 
of the lowest frequency phonon peak. It is found that the numerically measured renor- 
malization factor and sound velocity agree well with the theoretical predictions for FPU-/? 
chains, whereas obviously deviate from predictions for FPU-a/? chains. We stress that the 
asymmetric interactions play an important role responsible for breakdown of EPT, which 
provide stronger interactions among the phonons beyond the ones without momentum ex- 
change occurred in the chains with symmetric interactions. We also confirm that the analysis 
based on the lowest frequency phonon peak can obtain the renormalized phonon frequency 
in momentum conserved chain with asymmetric interactions where EPT break down. It is a 
challenge to correct theories to effectively renormalize anharmonicity induced by asymmetric 
interactions. 

Furthermore, the lowest frequency phonon peak analysis provide a test bed for the theories 
based on the long- wavelength approximations, such as the self-consistent mode coupling 
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theory 



lOj and hydrodynamics 12Hl4j. The current work together with our recent reports 



on heat conduction 



16 



17] have revealed the important roles of asymmetric interactions 
for phonon interactions in ID momentum conserved nonlinear chains, which have been 
underestimated in conventional theories. 
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